Which cosmological model — with dark energy or 
modified FRW dynamics? 



Marek Szydlowski 

Astronomical Observatory, Jagicllonian University, Orla 171, 30-244 Krakow, Poland 
Complex Systems Research Centre, Jagiellonian University, Reymonta 4, 30-059 
Krakow, Poland 

E-mail: szydloOoa . uj . edu . pi 

Wlodzimierz Godlowski 

Astronomical Observatory, Jagicllonian University, Orla 171, 30-244 Krakow, Poland 

Abstract. Recent measurements of distant type la supernovae (SNIa) as well as 
other observations indicate that our universe is in accelerating phase of expansion. 
In principle there are two alternative explanation for such an acceleration. While in 
the first approach an unknown form of energy violating the strong energy condition 
is postulated, in second one some modification of FRW dynamics is postulated. The 
both approaches are in well agreement with present day observations which is the 
manifestation of the degeneracy problem appearing in observational cosmology. We 
use the Akaike (AIC) and Bayesian (BIC) information criteria of model selection to 
overcome this degeneracy and to determine a model with such a set of parameters 
which gives the most preferred fit to the SNIa data. We consider five representative 
evolutional scenarios in each of groups. Among dark energy proposal the ACDM model, 
CDM model with phantom field, CDM model with topological defect, model with 
Chaplygin gas, and the model with the linear dynamical equation of state parameter. 
As an alternative prototype scenarios we consider: brane world Dvali Gabadadze 
Porrati scenario, brane models in Randall-Sundrum scenario, Cardassian models with 
dust matter and radiation, bouncing model with the cosmological constant and metric- 
affine gravity (MAG) inspired cosmological models. Applying the model selection 
criteria we show that both AIC and BIC indicates that additional contribution arises 
from nonstandard FRW dynamics are not necessary to explain SNIa. Adopting the 
model selection information criteria we show that the AIC indicates the flat phantom 
model while BIC indicates both flat phantom and flat ACDM models. 
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1. Introduction 

If we assume that the FRW model with a source in the form of perfect fluid well 
describes the present evolution of our Universe, then there is only one way of explain 
the observational fact that the Universe is accelerating [II El Ej -to postulate 
some gravity source of unknown nature which violates the strong energy condition 
Px + 3px > where px and px are energy density and pressure of dark energy, 
respectively [H E]. These different candidates for dark energy were confronted with 
observations [HllZlElinilinillllliailSllIlllISllinillZllIHl- While the most natural 
candidate for such a type of dark energy is vacuum energy (the cosmological constant A), 
we are still looking for other alternatives because the fine tuning problem with present 
value of the cosmological constant. In this context it is considered an idea of dynamical 
form of the equation of state (EOS) or decaying vacuum. As a first approximation one 

can consider the coefficient in the EOS wy = — in the form linearized around the 

A Px 

present epoch with respect to redshift z or the scale factor a: w(z) = u>o + W\z (or 
w(z) — Wo + wi(l — a), where a = 1 corresponds to the present epoch). The models 
with dust matter and such a form of the EOS for dark energy we called the dynamical 
EOS (DEOS) model. In the special case when W\ = and Wq < — 1 we obtain the CDM 
quintessence model or phantom CDM (PhCDM) model. Another interesting possibility 
of description dark energy offers a conception of description of dark energy in terms of 
Chaplygin gas 0. In this model the equation of state has a form p x = 

If we consider the FRW dynamics in which dark energy is present, the basic equation 
determining the evolution has the following form 

Peff k 
3 a 

where p e ff(a) is effective energy density of noninteracting "fluids", k = ±1,0 is the 
curvature index. Equation ((TJ can be presented in terms of density parameters 

H 2 

- m = n eff (z) + n kt0 (i + z) 2 (2) 

where ^ = j^, Q e ff(z) = Q m ^(l + z) 3 + Qx,of( z ) and Q m ,o is the density parameter 
for the (baryonic and dark) matter scaling like a -3 . For a = a (the present value of 
the scale factor) we obtain the constraint fi e //,o + ^fc,o — 1- 

We assumed that energy density satisfies the conservation condition 

Pi = -3H(pi +Pi), (3) 

for each component of the fluid p e ff = Ep^. Then from eq. (J2J) we obtain the constraint 
relation £jf2j j0 + ^fc,o — !■ 

All mentioned before directions which are coming toward to the description of dark 
energy in the framework of standard FRW cosmology can be treated as a representative 
approaches of explanation of the current Universe in terms of dark energy. In Table 
we complete all these models together with the dependence of Hubble's function H — - 
with respect to the redshift. We also denote the number of a model's free parameters 



# 2 =^-^. (1) 
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by d. Note that for the flat model, = 0, the number of the model parameters is 
equal d — 1. As the reference model we consider the flat FRW model with A = (the 
Einstein-de Sitter model with f2 mj0 — !)■ 

Since the discovery of acceleration of the Universe the theoretical and observational 
cosmology becomes in the state of permanent tension because of opposite aims of 
investigations. While the theoretical investigations go towards to generalization degree 
of consideration by adding some new model parameters, the observational cosmology 
tries to constraint these parameters. The main goal of observational cosmology is to find 
a model with a minimal number of essential model parameters. Then this, (of course 
it may be a naive model) is a starting point for the further analysis of constraints from 
the observational data. In the present observational cosmology such a role plays the 
concordance ACDM model and the PhCDM model [HI EDI ED E2j . 

Because nature of dark energy is unknown, it is considered another theoretical 
possibility that the phenomenon of accelerated expansion is actually a sign of a 
breakdown of the classical Friedmann equation which governs the expansion rate. In 
this context dark energy effects can be manifestation of a modification to the FRW 
equation arising from the new exotic physics. 

According to the brane world idea, the standard particles are confined on a hyper- 
surface which is called a brane, which is embedded in a higher dimensional bulk space 
time in which gravity could spread [23] . Then some additional contribution to the 
standard FRW equation may arise as a consequence of embedding our Universe in the 
higher dimensional bulk space. It is interesting that the cosmological models formulated 
in the framework of brane induced gravity can explain the acceleration of the Universe 
without the conception of dark energy. The additional term in the FRW equation drives 
the acceleration of the Universe at a late epoch when it is dominant [2H 1251 I2E] 1271 128j . 
Therefore the effect arising from existence of these additional dimensions can mimic 
dark energy through a modified FRW dynamics. As a prototype of evolutional scenarios 
arising from possible extra dimensions we consider two models: the Dvali Gabadadze 
Porrati (DGP) model and the Randall-Sundrum scenario of the brane-world (RSB) 
model. In the DGP model there is present a certain crossover scale r c that defines 
what kind of gravity an observer located on the brane observes. While for shorter 
then r c distances observer measures the standard gravitational force, for larger then 
r c distances gravitational force behaves like r~ 3 . In the RSB models there are present 
some additional parameters which are absent in the standard cosmology, namely brane 
tension A and dark radiation U (see Table [2J. If we assume the quantum nature of 
the Universe then quantum gravity corrections are important at both big bang and 
big rip singularities. Note that a big rip singularities can not be generic future state in 
phantom cosmology. Only if energy density is unbounded function of time one obtains a 
big rip/smash [29J. To account of quantum effects allows to avoid the initial singularity 
which phenomenologically can be modelled by a bounce [301 ED EU EH] El] • We consider 
bouncing models as a prototype of evolutional scenario in which quantum effect was 
important in its very early evolution. 
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Table 1. The five prototypes models explaining acceleration in terms of dark energy conceptions. 
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name of model 



H(z) 



free parameters 



1 DGP 
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Table 2. The five prototypes models explaining acceleration in terms of the modification of the FRW equation. 
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Freese and Lewis (35] have recently proposed to modify the FRW equation by 
adding a priori an additional term proportional to p™^ (dubbed the Cardassian term 
by the authors). If we consider a single fluid with energy density p then the Cardassian 
model stays equivalent to a two component noninteracting model. In the special case 
of dust matter there is a simple interpretation for the origin of this new "Cardassian 
term" — a perfect fluid satisfying the equation of state p = (n — l)p. Therefore, if n < 
the phantom cosmological models can be recovered. In our analysis we consider two 
fluids, matter and radiation. However, note that the radiation term is presently small 
with comparison to the matter term. We can assume = 0.0001 for radiation matter 



The next possibility is the MAG cosmological model based on a non-Riemannian 
gravity theory [HZ1ISH1- Because this model unifies both the RSB model with vanishing 
dark radiation and the model with spinning fluid 139 3 we include this model to the class 
of models which are based on the modified FRW equation. Because the RSB and MAG 
models without the cosmological constant do not fit well the SNIa data [101 EE1 HOj we 
also analyze the MAG and RSB models with the non- vanishing A term. The five models 
based on the modified FRW equation are presented in Table El 

Main goal of this paper is comparison of two conceptually different classes of 
models using the information criteria of models selection [IT]. The Bayesian information 
criterion gives important information whether additional parameters introduced by "new 
physics" are actually relevant and to have impact on the current Universe. 

2. Distant supernovae as cosmological probes dark energy or the modified 
FRW equation 

Distant type la supernovae surveys allowed us to find that the present Universe is 
accelerating PQI2|- Every year new SNIa enlarge available data sets with more distant 
objects and lower systematics errors. Riess et al. [Sj compiled the latest samples which 
become the standard data sets of SNIa. One of them, the restricted "Gold" sample of 
157 SNIa, is used in our analysis. 

For the distant SNIa one can directly observe their the apparent magnitude m 
and redshift z. Because the absolute magnitude M. of the supernovae is related to its 
absolute luminosity L, then the relation between the luminosity distance d^ and the 
observed magnitude m and the absolute magnitude M has the following form 



[36 



ra-M = 51og 10 d L + 25. 
Instead of dL, the dimensionless parameter Dl 



(4) 



D L = H d L 



(5) 



is usually used and then eq. (pQ) changes to 
pi = m - M = 5 log 10 D L + M 



(6) 
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where 

M = -5 log 10 H + 25. (7) 

We know the absolute magnitude of SNIa from the light curve. The luminosity 
distance of supernova can be obtain as the function of redshift 



where Qko = — -J? and 



(sinh(x) for < 
x for k = (9) 

sin(x) for fc > 

Finally it is possible to probe dark energy which constitutes the main contribution 
to the matter content. It is assumed that supernovae measurements come with the 
uncorrelated Gaussian errors and in this case the likelihood function L can be determined 
from the chi-square statistic £ oc exp(— x 2 /2) where 

I theor obs\2 

x 2 = E~ (10) 

i °i 

while the probability density function of cosmological parameters is derived from Bayes' 
theorem [T]. Therefore, we can perform the estimation of model parameters using the 
minimization procedure, based on the likelihood function. 

In the modern observational cosmology there is present so called the degeneracy 
problem: many models with dramatically different scenarios are in good agreement with 
the present data observations. Information criteria of the model selection [H] can be 
used to solve this degeneracy among some subclass of dark energy models. Among these 
criteria the Akaike information (AIC) [12] and the Bayesian information criteria (BIC) 

are most popular. From these criteria we can determine the number of the essential 
model parameters providing the preferred fit to the data. 

The AIC is defined in the following way 

AIC = -2\n£ + 2d (11) 

where £ is the maximum likelihood and d is a number of the model parameters. The 
best model with a parameter set providing the preferred fit to the data is that minimizes 
the AIC. 

The BIC introduced by Schwarz [33] is defined as 

BIC = -2\n£ + d\nN (12) 

where N is the number of data points used in the fit. While AIC tends to favour models 
with large number of parameters, the BIC the more strongly penalizes them, so BIC 
provides the useful approximation to full evidence in the case of no prior on the set of 
model parameters 
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The effectiveness of using these criteria in the current cosmological applications 
has been recently demonstrated by Liddle [33] who, taking CMB WMAP data [45J, 
found the number of essential cosmological parameters to be five. Moreover he obtained 
the important conclusion that spatially-flat models are statistically preferred to close 
models as it was indicated by the CMB WMAP analysis (their best-fit value is 
Q tott0 = £A,o = 1.02 ± 0.02 at la level). 

In the paper of Parkinson et al. 03] the usefulness of Bayesian model selection 
criteria in the context of testing for double inflation with WMAP was demonstrated. 
These criteria was also used recently by us to show that models with the big-bang 
scenario are rather prefers over the bouncing scenario |3E]- 

Please note that both information criteria values have no absolute sense and only 
the relative values between different models are physically interesting. For the BIC a 
difference of 2 is treated positive evidence (6 cLS cL strong evidence) against the 
model with larger value of BIC |37| l4*%j . Therefore one can order all models which 
belong to the ensemble of dark energy models following the AIC and BIC values. If we 
do not find any positive evidence from information criteria the models are treated as a 
identical and eventually additional parameters are treated as not significant. Therefore 
the information criteria offer the possibility of introducing relation of weak order in the 
considered class of analyzed models. 

The results of calculation of the AIC and BIC in the context of dark energy models 
are presented in TablesEKl I n Table|3]we show results for dark energy models considered 
for both flat and non flat cases without any assumed extra priors, while in Table 13] 
we presented results for models with the modified FRW equation. In general case 
the number of essential parameters in the cosmological models with dark energy is in 
principal two, i.e., H , TO| o- It means that the flat model is favored in the light of 
the information criteria. We can observe two rival models which minimize the AIC and 
BIC. They are the ACDM model and the phantom CDM (PhCDM) model. One can 
observe that both BIC and AIC values assume lower values for phantom models. It 
can be regarded as a positive evidence in favor of the PhCDM model. Basing on these 
simple and objective information criteria we obtain that SNIa data favored the models 
with the initial (big bang) and final (big rip) singularities. 

At first we analyze three flat models with two free parameters, i.e., the flat ACDM, 
TDCDM and PhCDM models. There is a significant difference between predictions of 
these models. The ACDM model prefers a universe with fi mj o close to 0.3, the PhCDM 
model favors a high density universe while the TDCDM model favors a low density 
universe. In Fig. and Fig. |2] we present values of the AIC and BIC for these and 
Cardassian models. If fi m> o < 0.22 then the information criteria favor the TDCDM 
model. For Q mfi G (0.22, 0.34), the ACDM is favored while for Q mfi > 0.34 the PhCDM 
model is preferred. The AIC allows also the Chaplygin gas model. However the BIC 
again prefers the ACDM model against the Chaplygin gas model. 

The similar analysis with the use of the information criteria is done in the case of the 
assumed prior f2 m = 0.3 49J (Table |3J). In this case the model with A is preferred over 
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the models with phantoms, that is in contrary to the results obtained in the previous 
case of no priors for Q m fl. It clearly shows that more precise measurements of O m) o will 
give us the possibility to discriminate between the ACDM and PhCDM models. 

For the model models with the modified FRW equation information criteria prefferes 
only the Cardassian model. However this model is preferred in the case of a high density 
universe, with f2 m0 > 0.44 which seems to be too high with comparison with present 
extragalactic data [49J. 

The comparison of PhCDM models with different fixed value of wx shows Fig. El 
and Fig. 0] Note that it is required to have higher matter content to equilize the greater 
negative effect of lower (more negative) value of the factor w. The model which minimize 
the values of both AIC and BIC is the model with w = —2.37. 

Fig. Eland Fig.Elshows the AIC and BIC in respect to Q m ,o, for low density universes 
(the TDCDM and DGP models) and the ACDM model. Dependent on value of fl m ,o the 
different cosmological models are selected. We find that TDCDM model is distinguished 
for Q mfi < 0.16, DGP model for fl mfi G (0.16,0.24) and ACDM for fi m>0 > 0.24. 

For deper analysis of statistical results it would be useful to consider the information 
entropy of the distribution which is defined as 

Entropy = -Ef { \og a (fi) (13) 

where a is a number of independent states of the system. 

In Table the value of entropy for four flat models (with topological defect, with 
the cosmological constant, with phantom and for the Cardassian model). The value of 
entropy for one dimensional PDF (fi m)0 ) is also presented. We can see that in both 
cases we obtain minimal value of entropy for the PhCDM model. 

Basing on these simple and objective information criteria and the minimum entropy 
principle we obtain that SNIa data favor the models with dark energy rather than based 
on modified FRW equations. Among the dark energy models the best candidates are the 
PhCDM and ACDM models. However, to make the final decision which model describes 
our Universe it is necessary to obtain the precise value of f2 m> o from independent 
observations. 

3. Conclusion 

The main goal of this paper is to decide which class of models: dark energy or models 
based on the modified FRW equation are distinguished by statistical analysis of SNIa 
data. For this aims we use the Akaike and Bayesian information criteria. 

It was considered two groups of five different models containing dust and dark 
energy (the first group) and dust matter and an additional term which modifies FRW 
dynamics (the second group). Our main conclusion is that both criteria weigh in favor 
of the flat dark energy models. We argue that there is no strong reason to favor models 
with modified FRW dynamics over FRW dark energy models. 

We are also able to decide which model with dark energy is distinguished by the 
statistical analysis of SNIa data. To do this we use the Akaike and Bayesian information 
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criteria. The former criterion weighs in favor of the flat phantom model while the latter 
distinguishes the flat phantom and ACDM models. Assuming the prior fi m = 0.3 both 
the AIC and BIC criteria weighs in favor of the model with dark energy, namely the flat 
ACDM model. 

The further conclusions are the following. 

• The minimal number of essential parameters in the cosmological models with dark 
energy is in principal two, i.e., (H ,ft mt0 ). The list of essential parameters may be 
longer, because some of them are not convincingly measured with present data, like 
the parameter W\. 

• The curvature density parameter does not belong to the class of essential parameters 
when all the rest parameters are without any priors (with no fixed fl m ,o)- At this 
point our result coincides with analogous result obtained by Liddle who found it 
basing on other observations. 

• If we consider models in which all model parameters are fitted then the PhCDM 
model with double initial and final singularities is distinguished. 

• When we consider the prior on f2 mj o then while for Q m ,o < 0.16 the model with 
two-dimensional topological defect is favored. The value fl m ,o €= (0.16; 0.24) favor 
DGP model while the value fi m>0 e (0.22; 0.34) favor the ACDM model. For 
^m,o G (0.22; 0.34) the phantom model is preferred. With high density Universe 
(f2 m > 0.44) Cardassian model (with modified FRW equation of state) is preferred. 

During our analisis we have use recently available SNIa data and of course future 
SNAP data will provide mor sophisticated information to distinguish between FRW 
dark energy models and the models with modified FRW dynamics. 

To make the ultimate decision which model describes our Universe it is necessary 
to obtain the precise value of fl m ,o from independent observations. 
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Table 3. The values of the AIC and BIC for the models from Table □ both for flat 
and non-flat cases. 
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Table 4. The values of the AIC and BIC for the models from Table H both for flat 



and non-flat cases. 
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Table 5. The values of the AIC and BIC for the models from Tabled with the prior 
^m,o = 0.3 both for flat and non-flat cases. 
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Figure 1. The values of the AIC in respect to fixed value fi m ,o for four flat models 
(with topological defect, with the cosmological constant, with phantom and for the 
Cardassian model). 



Table 6. The values of the AIC and BIC for models from Table [5] with the prior 



^m.o — 0.3 both for flat and non-flat cases. 
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Figure 2. The values of the BIC in respect to fixed value fl m ,o for four flat models 
(with topological defect, with the cosmological constant, with phantom and for the 
Cardassian model). 



Table 7. The value of information entropy for four flat models (with topological 
defect, the cosmological constant, phantom and for the Cardassian model). The value 
of entropy for the one dimensional PDF of f2 m .o is also presented. 
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